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1 Introduction and overview 

The standard model of particle physics is a local gauge-invariant quantum field theory in which 
symmetries play a fundamental role that includes the dependence of system properties under 
specific transformations such as charge conjugation (C), parity or space reflection (P) and time 
reversal (T). These individual symmetries and the combined CP symmetry are known to be 
violated in weak interactions, but the CPT combination appears to be conserved in nature U. 
A major consequence of CPT conservation is that the mass of any particle must equal that of 
its antiparticle. We focus on a measurement of the mass difference between the top and antitop 
quark. Since quarks carry color charge and hadronize into colorless particles before decaying, 
they cannot be observed as free quarks. The lone exception is the top quark, which due to 
its short lifetime decays before hadronization. The mass difference between the top quark 
and its antiquark was measured previously by the DO and CDF experiments, and showed no 
significant deviation from zero ||2]-2|. 

This letter reports a measurement of the difference between the mass of the top quark (t) and of 
its antiparticle (t), with significantly reduced uncertainties, using tt events produced in proton- 
proton collisions at ^/s = 7 TeV, recorded with the Compact Muon Solenoid (CMS) detector 
at the Large Hadron Collider (LHC) |5|. We select events where one W boson, either from the 
top or antitop quark, decays into qq' (t bW^ — > bqq', or its charge conjugate), and the 
other W decays leptonically (t — > bW^ — > hi^V£, or its charge conjugate), where the lepton i 
is a muon or an electron. The data are split into £^ and samples that contain, respectively, 
three-jet decays of the associated top or antitop quarks. For each event category, the Ideogram 
likelihood method |6| is used to measure the mass of the top quark (mt) or antitop quark (m^), 
and the difference between the masses in the two categories of lepton charge is taken as the 
mass difference Amt = mt — m^. The Ideogram method was used previously (ZllHl to measure 
the mass of the top quark. The procedure incorporates a kinematic fit of the events to a tt 
hypothesis that is modified specifically for this analysis to consider only the top or antitop 
quark that decays to three jets. 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam- 
eter, providing a magnetic field of 3.8 T. The field volume houses the silicon-pixel and silicon- 
strip trackers, a crystal electromagnetic calorimeter (ECAL) and a brass /scintillator hadron 
calorimeter. The inner tracker reconstructs charged-particle trajectories within the pseudora- 
pidity range | // 1 < 2.5, where the pseudorapidity is defined in terms of the polar angle 6 relative 
to the counterclockwise-rotating proton beam as rj = — In (tan 0/2). The tracker provides an 
impact parameter resolution of ~ 15 jim and a transverse momentum (pj) resolution of ~ 1.5% 
for 100 GeV particles. The energy resolution is < 3% for the electron energies in this analysis. 
Muons are measured for \tj\ < 2.4 using gaseous detection planes based on three technologies: 
drift tubes, cathode-strip and resistive-plate chambers. Matching outer muon trajectories to 
tracks measured in the silicon tracker provides a transverse momentum resolution of 1 — 6% 
for the Pj values in this analysis. In addition to the barrel and endcap detectors, CMS has ex- 
tensive forward calorimetry. A more detailed description of the CMS detector can be found in 
Ref. 0. 
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4 Event reconstruction and selection 



3 Data and simulation 

This analysis is based on a data sample corresponding to an integrated luminosity of 4.96 ± 
0.11 fb^^ collected in pp collisions at a center-of-mass energy of 7TeV and recorded with the 
CMS detector. Events are selected through a trigger requiring an isolated electron or muon 
with Pj > 25 or 17GeV, respectively accompanied by at least three jets of pj > 30 GeV in each 
event. The acquired data are compared to a set of simulated pp collisions at v^s = 7 TeV. Most 
signal and background events are generated with the matrix-element generator MadGraph 
4.4.12 flOl, interfaced to pythia 6.4.22 [11 J for the parton showering, where tt events are gen- 
erated accompanied by up to three extra partons. The mlm algorithm [12J is used for matching 
the matrix-element partons to their parton showers. Singly produced top-quark events are gen- 
erated with the POWHEG event generator Ill3| and generic multijet events with pythia. The 
simulation of multijet events is used just to normalize a multijet-enriched control sample of 
data needed in the analysis (described below). The simulation also includes effects of pileup in 
pp collisions, which refers to additional pp interactions that can occur during the same bunch 
crossing or in those immediately preceding or following the primary generated process. The 
simulated event samples are normalized to the theoretical cross section for each process, as 
calculated with FEWZ fT4'| for W and Z production, with pythia for multijet production, and 
MCFM 1 15 1 for all other contributing processes. The generated events are then passed through 
the full CMS detector simulation based on GEANT4 [16J, and eventually reconstructed using 
the same algorithms as used for data. 

4 Event reconstruction and selection 

Events are reconstructed using the particle-flow (PF) algorithm illTl , which combines the infor- 
mation from all CMS sub-detectors to identify and reconstruct individual particles produced 
in the proton-proton collision. The reconstructed particles include muons, electrons, photons, 
charged and neutral hadrons. Muons are reconstructed using the combined information from 
the silicon tracker and muon system HTSl . Electron reconstruction starts from energy deposi- 
tions in the ECAL, which are then matched to hits in the silicon tracker and used to initiate 
a track reconstruction algorithm. This algorithm takes into account the possibility of signif- 
icant energy loss of the electron through bremsstrahlung as it traverses the material of the 
tracker ||T9 | . Charged particles are required to originate from the primary collision vertex, iden- 
tified as the reconstructed vertex with the largest value of Epj for its associated tracks. The list 
of charged and neutral PF particles originating from the primary collision vertex is used as 
input for jet clustering based on the anti-A:T algorithm Il20ll with a distance parameter of 0.5. 
Particles identified as isolated muons and electrons are excluded from jet clustering. The mo- 
mentum of a jet is determined from the vector sum of all particle momenta in the jet, and from 
simulation found to be typically within 5 — 10% of the true jet momentum. Jet-energy-scale cor- 
rections are applied to all the jets in data and simulation. Jets in data have a residual correction 
that is determined from an assumed momentum balance in dijet and photon+jet events. These 
corrections are defined as a function of pj and t] of the reconstructed jet so as to obtain a more 
uniform energy response at the particle level, which tends to equalize the jet response in data 
and simulation [21J. The energy of jets is also corrected for the presence of additional pileup 
from neutrals, as the neutral component of pileup is still present after rejecting the contribution 
from charged hadrons. 

Events in the f/+jets channel are required to contain only one isolated muon with pj > 20 GeV 
and \t]\ < 2.1, while the e+jets channel requires only one isolated electron with pj > 30 GeV 
and < 2.5. The relative isolation is calculated from the other PF particles within a cone of 
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AR = Y^(A//)2 + (A^)2 < 0.4 around the axis of the lepton, with cp representing the azimuthal 
angle. It is defined as I^ei = (Wrged + ^photon + Jneutrai)/pT/ where /charged is the transverse 
energy deposited by charged hadrons, and /photon and /neutral are the respective transverse en- 
ergies of photons and neutral particles not identified as photons. Leptons are considered to 
be isolated when /j-ei < 0.125. Furthermore, events must have at least four reconstructed jets 
with Pj > 30GeV and < 2.4. Additional event selection criteria are discussed in Section |5] 
Table [T] gives the number of events observed in data following all selections, and the number 
expected from simulation, separately for events with fi^, e+ and e^. The kinematic char- 
acteristics of the multijet background are estimated from a data sample of events that pass all 
selections, but with an inverted lepton-isolation criterion of /^i > 0.2. The number of expected 
multijet events passing all selections is normalized to the PYTHIA simulation. 

Table 1: Number of events following full selection of /i++jets, //^+jets, e++jets and e~+jets 
events in data, and the expectation from simulations before any rescaling. Uncertainties are 
purely statistical and do not include contributions from production cross sections, integrated 
luminosity, detector acceptance, or selection efficiencies, as discussed in Refs. ||221|23| . 



Sample 


^++jets 


]d +jets 


e++jets 


e +jets 


tt 


15028 ± 56 


15006 ± 56 


10649 ± 47 


10611 ± 47 


W+jets 


11180 ± 149 


7314 ± 121 


7783 ± 125 


5523 ± 105 


Z/7*+jets 


1410 ± 25 


1516 ± 26 


1607 ± 27 


1685 ± 27 


Single top 


951 ± 7 


850 ±7 


675 ±6 


610 ±6 


Multijet 


483 ± 90 


196 ± 57 


722 ± 246 


1413 ± 485 


Total 


29050 ± 185 


24882 ± 147 


21436 ± 281 


19842 ± 499 


Observed 


27038 


23928 


22999 


21111 



Agreement between data and simulation in the number of selected events (normalization) is 
less important for this analysis than agreement for their kinematic distributions. The simulated 
signal and background events are therefore rescaled through a single global factor to match 
the number of events observed in data, keeping the relative background fractions fixed to the 
expectations from simulation. After this rescaling, a comparison of simulation and data for 
several key distributions is shown in Fig. [l] In general, the data appear to be well modeled 
by the simulation. The small possible deviations between data and simulation at large jet px 
values have little impact on this analysis, as most tt events have jet transverse momenta below 
200 GeV. 

5 Kinematic fit 

A kinematic fit of £+jets final states to a tt hypothesis, applying the constraints of transverse- 
momentum conservation, the assumed equality of nit and nii, and the accepted value of 80.4 GeV 
for the mass of the W boson (mw)/ has been one of the successful methods for extracting the 
mass of the top quark from tt events. The basic features of this type of kinematic fit are de- 
scribed in Refs. ||24ll25| . The fit we use corresponds to a modification that reconstructs the mass 
of the three-jet decays of top quarks (t — > bW — ?> bqq') by varying the momenta of the two 
jets that are assigned to the W boson, using niyj as a constraint. For each event, the four jets 
with highest transverse momentum (leading jets) are considered in the fit. These four jets can 
be associated with the four quarks for the tt-decay hypothesis tt — > bbW+W^ — > bbqq'fi/^ in 
24 possible ways. However, since the interchange of the two quarks from W-boson decay (qqO 
offers the same mass information, only 12 of these combinations provide unique solutions. The 
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5 Kinematic fit 
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Figure 1: Panels (a) and (b) display the transverse momenta of all jets in tt events, for and 
£^ events, respectively, while panels (c) and (d) give the number of selected jets per event 



that pass all selections. The simulation is normalized to the number of events observed in data. 
Overflows are included in the last bins of the distributions. The ratio of the number of observed 
events in data to the number of events expected from simulation is shown at the bottom of each 
plot. The error band corresponds to the systematic uncertainties related to jet energy scale, jet 
energy resolution, background estimation and modeling of pileup. 

four leading jets do not always originate from the quarks of the tt decay, because of the pres- 
ence of additional jets from gluon radiation. In simulated tt events, the three quarks arising 
from the top quark that decays to three jets are among the four leading jets in 70% of all such 
^+jets events. 

The kinematic fit is performed for each of the 12 jet-to-quark assignments. However, before 
implementing the fit, additional corrections are applied to correct jet energies to the parton 
level. These are derived separately for light-quark jets and for b jets in bins of and pj', 
by comparing the transverse energies Ej of selected jets with the Ej of generated partons in 
simulated tt events. The correction factors depend on the flavor of each jet for a given jet-to- 
quark assignment, and are about 4% larger for b jets than for light-quark jets. The parton- 
corrected jets used as input for the kinematic fit are parametrized by their Ej, 6, and (p. The 
resolutions of the reconstructed jet quantities are also used as input, and are obtained from the 
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width of the distributions for differences in Ej, 6, and (p between parton-corrected jets and the 
nascent parton values. As indicated above, the kinematic fit adjusts the momenta of the two 
light jets, taking their corresponding resolutions into account, while keeping the E/p of each 
jet fixed. Only solutions with x^/ '^dof < 10 ^re accepted, where n^of (= 1) corresponds to the 
number of degrees of freedom in the fit. An event is rejected if no combination of jets passes the 
X^/ ^dof requirement. The efficiency of this requirement in simulated tt events is 88%. The most 
important gain from the kinematic fit is that it improves the resolution on the top-quark mass. 
For correct jet combinations, the mass resolution is improved from ^ 15 GeV to ~ 10 GeV, as 
estimated from simulated tt events with nit = 172.5 GeV. 

The fitted values of the top-quark mass m,-, the uncertainty on the mass di and the xf, obtained 
for each combination of jets i, are used as input to the Ideogram method. A comparison of these 
variables between data and simulation is given in Fig. |2j for just the jet combination with the 
smallest in each event. Good agreement is observed between data and simulation. 

6 The Ideogram method 

In the Ideogram method, the mass of the top quark is measured using a likelihood defined as a 
function of Wt, as follows: 

Cevent{x;y \ nit) = /tt At(^;i/ I ^t) + (1 - ftt)Pbkg{x). (1) 

This equation expresses the likelihood for any mass nit in terms of a sum of probability densities 
from tt and background components. The fraction f^^ of the tt component is taken from Table [l] 
The functions Ptt(^' ¥ I ^t) and Pbkg(^) depend on the observables x, respectively for the tt and 
for the background hypotheses, where x includes the number of b-tagged jets n^,, the lepton 
charge q^, and the m, for each combination of jets i in the event. The quantities y represent the 
values of the parameters cr, and xf from the kinematic fits, and are used to parametrize P^t, as 
shown in Eq. (|4| below. The number of b-tagged jets is obtained using the Simple Secondary- 
Vertex High-Efficiency algorithm (SSVHE) Il26l at its "medium" working point. It is assumed 
that the background probability Pbkg can be described just by the probability density for the 
main background from W+jets. This is acceptable, as the contributions from other backgrounds 
are expected to be small, and their probability densities differ greatly from that from tt, but have 
distributions similar to that for W+jets events. 

Furthermore, it is assumed that the number of b-tagged jets and the lepton charge are uncorre- 
lated with the mass information in a given event. This means that the signal and background 
probabilities can be written as the product of a probability to observe n^, b-tagged jets, a proba- 
bility to observe a certain lepton charge q'', and a probability to observe Xmass/ which represents 
the set of mass variables m, in an event: 

Ptt{x;y I nit) = -Ptt("b) • Af('?0 • -Ptt(^mass;i/ 1 m); (2) 

^'bkg(^) = -Pbkg("b) • Phkgiq^) ■ -Pbkg(^mass)- (3) 

These probability densities for the number of b-tagged jets and lepton charge for signal, Ptt("b) 
and P^i{q^), and for background, Pbkg(wb) and Pbkgl'?^)/ are taken from simulation. The reason 
for including b-tagging at this point is to reduce the impact from non-tt background, while 
the reason for considering the probability distributions for lepton charge in background is to 
account for the dependence of W+jets and single-top events on the charge of the lepton. 

The tt probability for each event contains two terms, one representing the probability that a 
jet combination has the correct jet-to-quark assignment, and the second term expressing the 
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6 The Ideogram method 
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Figure 2: Panels (a) and (b) sh ow the dist ributions in fitted top-quark mass for the smallest fit- 
values, which are given in (c) and (d) for the kinematic fits for £^"+jets and i'^+jets events. 



respectively. The simulation is normalized to the number of events observed in data. The last 
bins include the contributions from overflow. The ratio of the number of observed events in 
data to the number of events expected from simulation is shown at the bottom. The error band 
corresponds to the systematic uncertainties related to jet energy scale, jet energy resolution, 
background estimation and modeling of pileup. 

probability that a jet combination has a wrong jet-to-quark assignment, which, summed over 
all possibilities i in each event, yields: 



-Ptt(^mass;j/ I mt) = J^^^i ^/gc J dm'G{mi \ m',ai)B{m' \ ni^Jt) + (1 -/gc)W(m, | mt)j . 



(4) 



The parameter /gc reflects the probability that the jet combination with highest weight Wj (de- 
fined below) corresponds to the correct jet-parton matching, as obtained from tt simulation, 
separately for events with Wb = 0, 1, and > 1. The probability for the correct jet combination 
is defined by the convolution in Eq. Q of a Gaussian resolution function G(m, | m',ai) and 
a relativistic Breit-Wigner distribution B{m' \ mt,Tt). The width of the top quark Ff is fixed 
to 2 GeV. The Gaussian function describes the mass resolution for each jet combination. It is 
centered at the Breit-Wigner-distributed value of the top-quark mass (m') and has a standard 
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deviation equal to the uncertainty on the fitted top-quark mass (cj-,). If the smallest xj in an 
event iXmm) larger than the value of n^of/ all the o-j for the event are scaled up by a factor 

\//Cmm/"dof- The symbol W(?«, | nit) in Eq. (4 1 represents the probabilities for the wrong jet 
combinations, which are parametrized using analytic functions fitted to the mass distribution 
of jet combinations from simulated tt events known to have wrong jet-to-quark assignments. 

The probability for tt signal of Eq. (jljl is calculated as a sum over all fitted jet combinations with 
X^ain < each weighted by: 

Wi = exp (^-^xfj^h- (5) 

The first factor above represents the likelihood for the kinematic fit with that combination of 
jets, while the second factor reflects the degree of compatibility with the observed b-tagging 
assignments: 

= n p'' 

where the index / runs over all jets considered in the fit, and the probabilities equal (1 — £/), 
£;,, or (1 — £(,), depending on the flavor assigned to each jet, and whether the jet is b-tagged. 
The b-tag efficiency (£f,) is 60.6 ± 2.5%, and is calculated from tt simulation using the scale 
factors between data and simulation and the corresponding uncertainties from Ref. 126]. The 
rate for tagging light-flavor jets (£/) is taken from Ref. [27J and equals 1.4 ± 0.3% for jets with 
50 < Pt < 80 GeV. The individual weights Wj are normalized to sum to unity for each event. 

The background probability Pbkg in Eq. ([Tjl does not depend on the mass of the top quark, 
and has only minimal dependence on the jet-quark assignments. The distribution is therefore 
defined by the mean of the combined distributions of all solutions for m, in simulated W+jets 
events, and fitted to an analytical function. 

The combined likelihood for the full event sample is calculated as the product of the individual 
event likelihoods for all selected events. The fitted top-quark mass and its statistical uncertainty 
are extracted from this combined likelihood. While /j^ can be treated as a free parameter of the 
fit 13, in this analysis it is fixed to the expected value (cf. Table [T| and the uncertainty on the 
signal fraction is taken into account as a systematic uncertainty of the method. 



7 Calibration of individual mass measurements 

The likelihood defined in Eq. ([T| for each event corresponds to a simplified model for the en- 
suing analysis, which means that the resulting combined likelihood reflects an approximate 
quantity. To correct for possible bias in the estimated mass or in the estimate of statistical 
uncertainty, a calibration of the procedure is performed using pseudo-experiments. In these 
pseudo-experiments, events are picked randomly from samples of simulated events represent- 
ing the major contributing processes in Table [ij implementing Poisson fluctuations around the 
respective means as expected in true data. The distributions for multijet events are modeled 
using control samples of data, as described in Sec. [l] For tt signal, nine samples of simulated 
events are generated for top-quark masses between 161.5 and 184.5 GeV. The calibration is 
performed for the accepted inclusive (> 3 jets) samples of £+jets events. 

The widths of pull distributions and the bias on the estimated top-quark mass as a function of 
generated mass are shown for the combined and events in Fig. |3] The pull is defined as 
the standard deviation of a Gaussian function fitted to the distribution of [nij — (m))/ crj, where 
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8 Measurement of the t — t mass difference 



nij is the estimated top-quark mass in each pseudo-experiment, aj its estimated statistical un- 
certainty and (m) the mean of the estimated top-quark masses over all pseudo-experiments at 
a fixed input mass. Since the standard deviation of the pull distribution appears to be slightly 
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Figure 3: (a) Width of the pull distribution and (b) bias on the extracted top-quark mass, as 
a function of generated top-quark mass for ^+jets events simulated in ensembles of pseudo- 
experiments. 

larger than unity, the statistical uncertainties on the final mass measurement are scaled up by 
that discrepancy (~ 16%). Also, as seen from Fig. |3| the bias on the estimated top-quark mass 
depends linearly on the generated top-quark mass. Although these biases are within 2 GeV for 
most of the range of interest, they are corrected using the fitted linear calibration given in the 
figure. The bias on the estimated top-quark mass and width of the pull as a function of gen- 
erated mass are shown separately for £++jets and £^+jets events in Fig. |4] (after implementing 
the inclusive ^+jets calibration from Fig. |3|. The results show that, within statistical precision, 
the separate £++jets and £^+jets events do not require additional independent corrections. 



8 Measurement of the t — t mass difference 

For the final measurement of the mass difference, we apply the analysis separately to ^^+jets 
events and to ^++jets events, and take the difference of the two extracted values. In the f/+jets 
channel, the individual measurements yield a mass difference of: 

Amt = 0.13 ± 0.61 (stat.) GeV, 

and in the e+jets charmel: 

Amt = -1.28 ±0.70 (stat.) GeV, 
and when the method is applied to the combined e+jets and f/+jets samples: 

Amt = -0.44 ±0.46 (stat.) GeV. 

The results for Amt are compatible with the expectation from the hypothesis of CPT symmetry, 
which forbids a mass difference between the top quark and the antitop quark. Also, the average 
fitted top-quark mass is found to be mt = 173.36 ± 0.23 (stat.) GeV, which is in agreement with 
previous measurements of mt Il28] - |3TI , even ignoring systematic uncertainties. 
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Figure 4 



9 Systematic uncertainties 



Many of the systematic uncertainties relevant for the absolute measurement of nit, such as 
the calibration of the overall jet energy scale, are reduced in the context of this measurement, 
as such systematic effects tend to alter the measured properties of top and antitop quarks in a 
similar and correlated manner. Several other sources of systematic uncertainty on the modeling 
of physical processes evaluated in absolute nit measurements are not expected to affect the 
measurement of Anit. These include modeling of hadronization, the underlying event, initial 
and final-state radiation, changes in factorization and renormalization scales, and the matching 
of partons to parton showers, and are not considered further in this analysis. 

Systematic uncertainties for other effects considered in the measurement of nit are included 
together with additional sources potentially relevant for Anit, such as lepton-charge identifica- 
tion and a possible difference in jet response to b and b quarks. These are listed in Table [2} and 
described in greater detail below. In all cases, the effects are evaluated using simulated event 
samples, by comparing the nominal sample to one where the systematic effect is varied by ±1 
standard deviation. Statistical uncertainties on the observed mass shifts are evaluated using 
the resampling technique of Ref. |32|, and are listed in Table |2j For most systematic uncertain- 
ties, the statistical significance of the observed shift in Anit is small. We therefore quote the 
observed shift as a systematic uncertainty when it is larger than the statistical uncertainty, and 
otherwise we quote just the statistical uncertainty. The total systematic uncertainty is taken to 
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9 Systematic uncertainties 



be the quadratic sum of the values quoted for each source. 

Table 2: Overview of systematic uncertainties on Amt. The total is defined by adding in quadra- 
ture the contributions from all sources, by choosing for each the larger of the estimated shift or 
its statistical uncertainty, as indicated by the bold script. 

Source Estimated effect (GeV) 



Jet energy scale 0.04 ± 0.08 

Jet energy resolution 0.04 ± 0.06 

b vs. b jet response 0.10 ±0.10 

Signal fraction 0.02 ± 0.01 

Difference in W+/W-production 0.014 ± 0.002 

Background composition 0.09 ± 0.07 

Pileup 0.10 ± 0.05 

b-tagging efficiency 0.03 ± 0.02 

b vs. b tagging efficiency 0.08 ± 0.03 

Method calibration 0.11 ± 0.14 

Parton distribution functions 0.088 

Total 027 



Overall jet energy scale. The uncertainty related to the overall jet energy scale is estimated by 
changing the energy of all jets within their pj and //-dependent uncertainties. This un- 
certainty contains contributions from the uncertainty on pileup and flavor dependence 
of jets. These are measured using 2010 data II2TII , and have several additional uncertain- 
ties added in quadrature: 1.5% from a mismatch in absolute energy between Z+jets and 
7+jets events, and 1.5% for jets with \r]\ > 1.3 from an //-dependence on the relative scale. 
Because top and antitop quarks at the LHC are produced with slightly different distri- 
butions in rapidity, an //-dependence for jet response can lead to a small residual effect 
on A/Wf While the average extracted top-quark mass shifts by as much as ±2.3 GeV, the 
observed effect on Amt is only 0.04 ± 0.08 (stat.) GeV. We quote the statistical precision on 
the shift as a systematic uncertainty of the measurement. 

Jet energy resolution. Previous measurements of jet energy resolution in data have indicated 
that it is 10% worse than in the simulation [21J. The resolution in simulated events used 
for calibration is therefore degraded accordingly. The uncertainty on this 10% depends 
on f], and equals ±10% for jets within \r]\ < 1.5, ±15% for jets within 1.5 < |//| < 2.0, 
and ±20% for jets with \t]\ > 2.0. Based on generated parton energies, the resolution for 
each jet is scaled up and down within these uncertainties. Half of the difference between 
such up and down changes yields a —0.04 ± 0.06 (stat.) GeV difference in A/t/f While this 
is expected to be the same for the nit arid measurements, a residual effect is possible 
through the asymmetry in the composition of the background. This possibility is included 
as a systematic uncertainty on Amt. 

Jet energy scale for b and b. A dedicated study is performed to assess the jet response for b 
and b jets, by comparing the reconstructed jet pj to the original parton pj in tt simu- 
lation as a function of jet // and pj. The pythia simulation describes differences in the 
fragmentation of b and b jets, including B-B oscillations, and the CMS detector simulation 
includes differences in calorimeter response for K+ and particles. As the PF algorithm 
reconstructs charged hadrons using tracks when available, the impact of such differences 
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in calorimeter response on the reconstructed jet energy is expected to be small. On av- 
erage, the ratio of b to b response is found to be 0.999 ± 0.001, compatible with unity. In 
principle, the calibration would correct for such differences, albeit with limited statistical 
precision (see below). Nevertheless, we quote 100% of the corresponding shift of 0.10 GeV 
as a systematic uncertainty on A/Wf 

Signal fraction. An incorrect fraction of tt signal events in the simulation would, in princi- 
ple, bias the calibration procedure. Changing the relative signal fraction by ±20%, while 
keeping the background composition fixed, yields an effect of ^ 0.02 ± 0.01 GeV on ArWf 

Difference in W^fW^ production. The difference in production cross sections of W+ and 

bosons in pp collisions leads to different levels of W+jets background and different back- 
ground composition in £++jets and £^+jets channels. This can affect the calibration pro- 
cedure and lead to a small bias in Anit. The measured inclusive W+/W^ ratio is in agree- 
ment with theoretical prediction within a precision of 3.5% II33II , and has been studied by 
CMS as a function of pseudorapidity |34|. Varying the and backgrounds by 2% 
in opposite directions, thereby affecting the relative ratio of W+ and by 4%, changes 
Anit by — 14 ± 2 MeV, which is quoted as the systematic uncertainty resulting from the 
difference due to unequal yields of W+ and . 

Background composition. To evaluate any residual effects related to distributions and com- 
position of the background, we investigate the effect of removing completely each source 
of background from the calibration procedure, while keeping the signal fraction con- 
stant. We quote 30% of the total shift observed in Amt when we remove W+jets (—0.26 ± 
0.20 GeV), Z+jets (0.05 ± 0.04 GeV) and single top-quark production (0.05 ± 0.02 GeV), and 
100% for the background from multijet events (—2 ± 5 MeV). For each contribution, we 
take the larger of the observed shift in Amt or its statistical uncertainty and add the four 
sources in quadrature. 

Pileup. The simulated events used in this analysis contain contributions from pileup, and are 
reweighted to match the estimated dependence of pileup on instantaneous luminosity. 
The systematic uncertainty is estimated by changing the mean value of the number of 
interactions by ±0.6, and taking the average of the two shifts in Amt as the systematic un- 
certainty. This covers the uncertainty in the modeling of pileup as well as the uncertainty 
on the calculation of event weights. The uncertainty on the weights is dominated by un- 
certainties on the total inelastic cross section and on the measured luminosity, both of 
which are used in the reweighting. To further investigate any additional effects related to 
high pileup conditions during high-luminosity running, the measured values are exam- 
ined as a function of the number of reconstructed primary vertices in Fig. |5] No adverse 
effects are observed, and the results for Amt are stable and statistically compatible with 
no dependence on the number of pileup events in the data. 

B-tagging efficiency. A shift in b-tagging efficiency can affect the impact of background pro- 
cesses on the Ideogram method. In particular, a difference in distributions of positively 
and negatively charged particles in the background can affect Amt. As indicated previ- 
ously, the SSVHE tagger is used primarily at its "medium" working point. To quantify 
the impact of b-tagging efficiency, we vary the threshold defining the working point, 
thereby producing a relative change in efficiency of ±4% p6) . The changes are applied 
in the same direction or opposite direction for the £''"+jets and £^+jets samples, and the 
corresponding shifts on Amt of 0.03 GeV (b-tagging efficiency) and 0.08 GeV (b versus b 
tagging efficiency) are quoted as systematic uncertainties. 
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9 Systematic uncertainties 
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Figure 5: Difference between (a) the measured in each bin and the average m^, and (b) the 
value of Amt both for data and for simulation, as a function of the number of reconstructed 
collision vertices. The results of fitting the data with a linear function are also shown, (c) 
Distribution for the number of reconstructed collision vertices in data and in simulation, after 
reweighting for pileup. 

Misassignment of lepton-charge. The leptons are used only in triggering and splitting the 
data into £++jets and £^+jets events, but not in mass reconstruction. A wrong assignment 
of charge can affect the calibration of Awt ^ a way that is not recovered in the overall pro- 
cedure. It can also lead to a cross-contamination of the two event samples, which can bias 
or dilute the measurement. The rate of charge misassignment in muon reconstruction is 
measured with cosmic muons [35] and collision data [34J to be of the order of 10^"' to 
10^^% in the transverse momentum range of this measurement. For electrons, the rate 
of charge misidentification ranges from 0.1% to 0.4%, depending on pseudorapidity ||34| . 
This means that the systematic uncertainty from charge misassignment is below 1% of 
the measured Amt value, which is negligible and is therefore ignored. 

Trigger. The trigger requires the presence of a lepton and at least three jets. As the lepton is 
not used in mass reconstruction, no systematic effect is expected from any mismodeling 
of the lepton trigger efficiency or px threshold. The requirement of three jets in the trig- 
ger is highly efficient for events with 4 jets with pj > 30 GeV. Any effect on kinematic 
distributions of the jets in selected events is therefore estimated to be small, and expected 
to affect the rrit and measurements equally. No uncertainty is quoted therefore for this 
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source. 

Method calibration. The effect is evaluated for simulated tt events at a mass of 172.5 GeV, 
showing a difference in mass bias between £++jets and £^+jets of —0.11 ± 0.14 GeV, which 
is statistically compatible with no effect. This confirms our expectation that there is no 
known effect in simulation that would lead to a difference in mass calibration between 
the two channels. Based on this observation, the combined £+jets calibration is applied 
both in the £++jets and the £^+jets channel. The statistical uncertainty on the calibration 
of the mass difference is quoted as a systematic uncertainty of 0.14 GeV. 

Parton distribution functions. The choice of parton distribution functions (PDF) can affect 
Amt, as they determine, for example, the difference in production of W+and W^, which 
is the dominant source of background. The simulated samples are generated using the 
CTEQ 6.6 PDF [36 1, for which the uncertainties can be described by 22 independent pa- 
rameters. Up and down changes in these parameters result in 22 accompanying PDF 
possibilities. Using a simulated sample of tt and background events, reweighted accord- 
ing to the deviation of each PDF from its original form, the sum of the larger shift ("up" or 
"down") for each change in PDF is taken in quadrature, to define an estimated combined 
uncertainty on A/«t of 88 MeV. 

10 Summary 

The mass difference between the top quark and the antitop quark, Amt = wJt — m^, is measured 
with the Ideogram method using the £+jets tt event sample collected by the CMS experiment, 
corresponding to an integrated luminosity of 4.96 ± 0.11 fb^^. This yields the result: 

Amt = -0.44 ± 0.46 (stat.) ± 0.27 (syst.)GeV 

The measured value is in agreement with the consequence of CPT invariance, which requires 
no mass difference between the top and antitop quarks. This is more precise by at least a factor 
three than any of the previous measurements. 
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